Spin-charge decoupling and the photoemission line-shape in one dimensional insulators 
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The recent advances in angle resolved photoemission techniques allowed the unambiguous experi- 
mental confirmation of spin charge decoupling in quasi one dimensional (ID) Mott insulators. This 
opportunity stimulates a quantitative analysis of the spectral function A(k,ui) of prototypical one 
dimensional correlated models. Here we combine Bethe Ansatz results, Lanczos diagonalizations 
and field theoretical approaches to obtain A(k,u>) for the ID Hubbard model as a function of the 
interaction strength. By introducing a single spinon approximation, an analytic expression is ob- 
tained, which shows the location of the singularities and allows, when supplemented by numerical 
calculations, to obtain an accurate estimate of the spectral weight distribution in the (k, uj) plane. 
Several experimental puzzles on the observed intensities and line-shapes in quasi ID compounds, 
like SrCuC>2, find a natural explanation in this theoretical framework. 

PACS numbers: 71.10.Fd, 79.60.-i 
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I. INTRODUCTION 



Since the theoretical prediction of the decoupling of 
spin and charge excitations in one dimensional (ID) 
models 1 , many experiments have long sought to verify 
this effect 2 . According to the spin-charge separation sce- 
nario, the vacancy (e + ) created by removing an electron 
in a photoemission experiment decays into two collective 
excitations (or quasi-particles) , known as spinon (s) and 
holon (h), carrying spin and charge degrees of freedom 
respectively. The recent observation of a well defined 
two-peak structure in the angle-resolved photoemission 
spectra (ARPES) of the quasi-lD materials SrCuC>2 and 
Sr2CuC>3 3 ' is deemed a significant clue of spin-charge 
decoupling, confirming previous expectations. 

However, other quasi one dimensional materials 5 fail 
to show distinct holon and spinon peaks, casting some 
doubt on the interpretation of ARPES experiments based 
on spin charge decoupling. A number of puzzling features 
also suggest that more physics, beyond the simple decay 
e + — ► s + h, is involved in the photoemission process: 
the spectral functions of SrCuC>2 and Sr 2 Cu03 reported 
by Kim et al. 3 and by Kidd et al. A systematically display 
broad line-shapes in contrast to the sharp edges expected 
on the basis of the available calculations on model sys- 
tems. The spectral intensity also appears considerably 
weaker in a half of the Brillouin zone, a feature often 
ascribed to cross section effects 6 . 

A quantitative theoretical understanding of ARPES 
in low dimensional systems is important and deserves 
a careful investigation because ARPES provides a direct 
experimental probe to the single particle excitation spec- 
trum, allowing for reliable estimates of the key param- 
eters governing the physics of strongly correlated elec- 
trons: the electron bandwidth and the Coulomb repul- 
sion. Here we will focus on the ID Hubbard model, a 
simple lattice model defined by just two coupling con- 
stants: the nearest neighbor hopping integral t and the 



on-site Coulomb repulsion U: 
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h.c. 



+ U^2n t1 n lL (1) 



Although several other terms, such as next-nearest hop- 
ping, further orbital degrees of freedom, temperature, 
disorder or lattice instabilities, would be necessary in a 
realistic model of these materials, we believe that an ac- 
curate investigation of the simplest hamiltonians should 
be performed before facing more challenging problems. 

The theoretical studies aimed at the investigation of 
the spectral properties of one dimensional models are 
either fully numerical, like Lanczos diagonalizations 2 
and density matrix renormalization group (DMRG) 
techniques 7 , or are carried out in the limiting cases 
of infinite 8 or vanishing 9 interaction U/t. In the for- 
mer case, they suffer from severe finite size effects, in 
the latter the interplay between charge fluctuations and 
strong correlations is not satisfactorily taken into ac- 
count. Monte Carlo studies of dynamical properties of 
quantum systems are instead hampered by the necessity 
to perform an analytic continuation to real times. 

In this paper, we provide the first quantitative evalua- 
tion of the full spectral function A(k,ui) of the ID Hub- 
bard model at half filling for intermediate and strong 
coupling U/t w . A formalism based on the Bethe Ansatz 
solution 11 , and supplemented by Lanczos diagonaliza- 
tions, is developed and is shown to provide a transparent 
description of the dynamical properties of mobile charges 
in Mott insulators. From this analysis we find that the 
ID Hubbard model does indeed contain the physics re- 
quired for a quantitative interpretation of photoemission 
experiments. In particular: i) the underlying free elec- 
tron Fermi surface plays a key role in defining the shape 
and the intensity of the ARPES signal, up to fairly large 
effective couplings U/t; ii) the power-law singularities 
which characterize the spectral function in one dimen- 
sion give rise to intrinsically broad peaks, whose width is 
proportional to the intensity of the line; Hi) ARPES data 



2 



are extremely sensitive to the Hubbard parameters and 
allow for a direct determination of the effective coupling 
constants in quasi ID materials. As a working example, 
we apply our method to SrCu02, where accurate ARPES 
data are available 2 , and we derive reliable estimates for 
t and U, 

The plan of the paper is as follows. In Section II we 
introduce and motivate the single spinon approximation 
which lies at the basis of our method, deriving the pre- 
dicted formal structure of the spectral function in one 
dimensional models. Section III shows how Lanczos di- 
agonalizations provide a precise quantitative estimate of 
the quasi-particle weight required for the evaluation of 
the spectral function. Then, in Section IV we discuss 
the weak coupling limit, where a thorough field theoret- 
ical analysis is available. The application to the case of 
SrCuC>2 is performed in Section V, while in the Conclu- 
sions we briefly discuss the generalization of our method 
to more complex one dimensional hamiltonians. 

II. THE ANALYTICAL STRUCTURE OF THE 
SPECTRAL FUNCTION 

The dynamical properties of one (spin down) hole in 
the half filled Hubbard model are embodied in the spec- 
tral function A(k,w) which, at zero temperature, can be 
written as: 

A(k,u,)= \(^i\c k ,i\^o)\ 2 S(uJ-E 1 +E ) (2) 
{|*i>} 

where |^o) is the ground state of the model at zero dop- 
ing, i.e. when the number of electrons N equals the 
number of sites of the lattice L, Eq is the correspond- 
ing energy and {| v I'i)} represents a complete set of one- 
hole intermediate states, of energies E±. The whole en- 
ergy spectrum of the Hubbard hamiltonian (1) can be 
obtained from the Lieb and Wu equations 11 . In the ther- 
modynamic limit, its structure has been thoroughly in- 
vestigated in a series of papers by Woynarovich 12,13 (see 
also the comprehensive book by Essler et al. Ref. 14 ). 
In summary, the exact excitation spectrum at half fill- 
ing and at arbitrary coupling U/t depends on two sets 
of "rapidities" describing the charge and spin degrees 
of freedom respectively. The excitation energy is al- 
ways written as the sum of contributions involving just 
two elementary excitations, representing collective quasi- 
particlcs: "holons" (of momentum kh and energy e/ l (fc? l )) 
and "spinons"(of momentum Q € (■!>%") an d energy 
e s (Q)). The simplest physical excitation created by the 
removal of an electron of momentum k gives rise to one 
holon and one spinon satisfying the momentum conserva- 
tion equation k = kh+Q- The total energy of this state is 
Ei = E +eh(kh)+e s (Q)- Besides this suggestive "decay" 
mechanism of the electron, other excited states also ap- 
pear in the exact spectrum: they are either multi-spinon 
and multi-holon states, or states involving the creation 
of double occupancies 13 . However, it is remarkable that 



the full excitation spectrum can be always expressed in 
terms of €h{kh) and e s (Q), showing that spin charge de- 
coupling holds, in the Hubbard model, at all values of 
U/t and at all energy scales 14 . The two quasi-particles, 
holon and spinon, are both collective excitations involv- 
ing an extensive number of degrees of freedom and can 
be approximately related to simple real space pictures of 
a "hole" and an unpaired spin only in the strong coupling 
limit, where spin-charge decoupling acquires a more in- 
tuitive meaning. As U —> the holon and spinon bands 
reduce to simple analytical forms 11 , closely related to the 
free particle band structure: e/ l (fc/ 1 ) = At cos(fch/2) and 
e s (Q) = 2t|cosQ|. 

While the whole energy spectrum of the Hubbard 
hamiltonian is known in detail, the matrix elements ap- 
pearing in Eq. (2) are of difficult evaluation. Moreover 
the summation over the intermediate states formally in- 
volves a number of terms exponentially large in N, mak- 
ing the exact implementation of the definition (2) im- 
practical. Our approach, which allows for the evaluation 
of the full spectral function in the thermodynamic limit, 
is based on the single spinon approximation: i.e. we ne- 
glect the contribution to the spectral function coming 
from all multi-spinon excited states and all excitations 
with complex rapidities, but we evaluate exactly the ma- 
trix elements involving one holon and one spinon. The 
accuracy of this method is tested a posteriori by use of a 
completeness sum rule and can be estimated of the order 
of few percents. Such a remarkable performance of the 
single spinon approximations is not unusual in one di- 
mensional physics: a known example is provided by the 
Haldane-Shastry spin model (HSM) 15 , where each inter- 
mediate state contributing to the dynamical spin corre- 
lation is completely expressible in terms of eigenstatcs 
of the HSM with only two spinons. In this case, only a 
small O(L) number of cigenstates contribute to the exact 
dynamical spin correlation function as proved in Ref. 16 . 
Similarly, in our approach, the most relevant intermedi- 
ate states are expressible is terms of cigenstates of the 
Hubbard model with only one spinon and one holon ex- 
citations. 

A first clue on the structure of the spectral function 
in one dimensional models can be obtained by analyzing 
the U —> oo limit, where double occupancies are inhib- 
ited and several exact results are available 8 . At half fill- 
ing (TV = L) the Hubbard hamiltonian is mapped onto 
a Heisenberg hamiltonian: each site is singly occupied 
and the ground state is a non-degenerate singlet of zero 
momentum 17 . When a hole of momentum kh is created, 
all the eigenfunctions of the Hubbard hamiltonian (with 
periodic boundary conditions) can be written as 18 : 

l*i >=-j= £ e ik ^4> H { yi ,...,y M )\x,{ yi }) (3) 

where \x, {yi}) represents the configuration of L — 1 elec- 
trons defined by the positions of the M = L/2 spin 
up ({yi}) and of the hole (x). The amplitude </># is a 
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generic eigcnfvmction of the Hciscnberg hamiltonian on 
the "squeezed chain" , i.e. on the L — 1 site ring defined 
by all the sites occupied by an electron. The intermedi- 
ate states |^i > entering the spectral function (2) have 
momentum — k relative to the ground state at half fill- 
ing. Due to the factorizcd form of the eigenfunctions (3) 
the total momentum of the state satisfies k = kh + Q 
where Q is the momentum of the Heisenberg eigcnfunc- 
tion (f>H, expressed in integer multiples of 2ir/(L — 1) 
in finite chains. In the thermodynamic limit the energy 
of the intermediate state is E\ = E Q + e/ l (fc/ l ) + e s (Q), 
where, to lowest order in J = 4t 2 /U, the first (holon) 
contribution is just the kinetic energy of a free particle 
(eh(kh) = 2t cos kh) and the second one (spinon) is the 
energy of the cigenstatc 4>h referred to the ground state 
energy of the Heisenberg ring of L sites 19 . This analysis 
shows, in an intuitive way, the origin of momentum and 
energy conservation in the decay process of the vacancy 
and suggests that, in the U —> oo limit, the most rele- 
vant contributions to the sum of intermediate states in 
Eq. (2) come from the lowest energy eigenstates <f>H of 
the Heisenberg hamiltonian for the L — 1 allowed mo- 
menta Q = -j^n (n = . . .X — 2). Accordingly, the sum 
over an exponentially large set of eigenstates {| , 3 / i >} in 
Eq. (2) can be (approximately) replaced by a sum over 
L — 1 single spinon states. This special set of intermediate 
states l^i) , which we argue provides the dominant contri- 
bution to the spectral weight for each spinon momentum 
Q, will be referred to as \ — k,Q) in order to emphasize 
the two quantum numbers which uniquely identify them. 
The single spinon approximation can be easily tested in 
the U — > co limit 8 where it proves extremely accurate. In 
the next Section we will show that it remains fully satis- 
factory also at finite coupling. In fact, it is known 14 that 
the eigenstate structure of the Hubbard model displays 
a remarkable continuity in U/t, the only singular point 
being the (trivial) free particle limit [7 = 0. However, 
when charge fluctuations are allowed for, by lowering the 
strength of the on-site repulsion U, the identification of 
the single spinon states \ — k,Q) is not easy, because the 
spinon momentum Q is not a good quantum number any 
more, although it can be still formally defined on the ba- 
sis of the Bethe Ansatz solution of the Hubbard model 13 . 
The key observation, which will be exploited in the next 
Section, is that the correct single spinon states can be 
identified at finite U via Lanczos or DMRG calculations 
by following adiabatically the evolution of the Heisenberg 
states as U is gradually decreased. 

Keeping only the single spinon states in the summation 
of Eq. (2), the spectral function in the thermodynamic 
limit becomes 

A(k,u) = f ^Z k {Q)5{Lo-e h {k-Q)-e s {Q)) 



matrix clement 



-E 



t ^ \v h (Q* -k)+v s (Q*)\ 



(4) 



Z k (Q)= lim (i-l)|(-fc,Q| Cfcji |v& )| 5 



(5) 



The sum in Eq. (4) runs over all the solutions Q*(k,ui) 
of the algebraic equation 



w = e h (k - Q) + e s (Q) 



(6) 



where th(kh) [ e s(Q)] is the known holon [spinon] exci- 
tation energy 13 and v h {k h ) = [v s (Q) = fy] the 
associated velocity. Equation (4) is the main result of 
this work: an explicit and computable expression for the 
spectral function of one dimensional models. In the spe- 
cial case of the Hubbard model, the Bethe Ansatz so- 
lution directly provides spinon and holon dispersions in 
the thermodynamic limit further simplifying the evalu- 
ation of the spectral function. Due to the presence of 
a spinon Fermi surface, the dispersion relation e s (Q) is 



defined only in the interval ^ < 



Where we have defined the quasi-particle weight as the 



< ^, 20 it van- 
ishes at the boundaries and has a single maximum at 
Q — 7r, while £h(kh) is an even and periodic function 
in the whole range — tt < kh < tt with maximum at 
k = 12,13 . The only missing ingredient in Eq. (4) is 
the quasi-particle weight Z k {Q) which defines the line- 
shape and intensity of the spectral function. Previous 
studies 21 have shown that in spin isotropic models, like 
the Hubbard model, the quasi-particle weight is a regu- 
lar function with square root singularities at the spinon 
Fermi surface Q = 7r±7r/2. This implies that A(k,ui) has 
power law singularities too, whenever cither Q* defined 
by Eq. (6) lies at the spinon Fermi surface, or when the 
total excitation velocity Vh(Q* — k) + v s (Q*) vanishes. In 
both instances, square root divergences are expected 22 : 
in the former case the location of the singularity iden- 
tifies the holon dispersion via (6) u> = €h(k + tt ± J); 
in the latter case the singularity is trivially due to band 
structure effects and does not necessarily corresponds to 
a pure spinon contribution as often assumed. However, 
at small to moderate interactions U/t, the holon veloc 
ity |w/i(fch)| displays an abrupt drop around kh ~ tt 
placing the band lower edge close to Q ~ tt + k, i.e. at 
w ~ £h(7r) + € s (tt + k), thereby following the spinon band 
for < k < T j. This particular feature of the Hubbard 
model dispersion is apparent in the shape of the holon 
spectrum 12 which sharply bends at kh ~ ±7r so to dis- 
play a vanishing charge velocity at band edges. This also 
agrees with the "relativistic" form of the holon spectrum 
predicted by bosonization at weak coupling 9 , as reported 
in Eq. (8). The expected location of the square root sin- 
gularities of the spectral function in the (fc, u>) plane is 
shown for few values of the coupling in Fig. 1. The holon 
branch (shown as full circles in the figure) marks precisely 
the holon excitation spectrum eh(kh) while the location 
of the singularities due to the band structure (shown as 
crosses in the figure) differs from the spinon e s (<5) dis- 
persion by less than O.lt. Note also that the curvature 
of the "spinon branch" displays a significant dependence 
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FIG. 1: Location of the singularities of the spectral function 
in the thermodynamic limit for three values of U in the plane 
(k,E = — uj). The (green) circles correspond to the singulari- 
ties of the Zk(Q) (holon branch), while the black stars to the 
extrema of the excitation spectrum. 




FIG. 2: Panel (a): holon bandwidth W h = e h (ji/2) - e h (-K) 
as a function of U/t. Panel (b): ratio between the spinon 
bandwidth W s = £s{k) and Wh as a function of U/t. 



on U/t, allowing for a rather precise experimental deter- 
mination of the effective coupling ratio. Therefore we 
conclude that precise photoemission data, able to iden- 
tify the singularities of the spectral function, do provide 
direct information on both holon and, within a good ap- 
proximation, also spinon excitations. 

The full holon bandwidth is always 4i at all couplings, 
due to the particle-hole symmetry of the Hubbard model 
but the upper and lower branches of the holon band are 
not symmetrical at finite U . This observation is relevant 
for the correct interpretation of photoemission experi- 
ments, because an estimate of the effective hopping inte- 
gral t is usually performed by measuring the half band- 
width of the upper holon branch 23 leading to a sizable 
overestimate of t. In Fig. 2 we show the bandwidth Wh 
of the upper holon branch (i.e. e^(7r/2) — e^(7r)) and 
the ratio between the spinon and the holon bandwidths 
W s /Wh as a function of the coupling U/t. Both quanti- 
ties, which allow for a direct estimate of t and U/t from 
ARPES, show a remarkable (even non monotonic) depen- 
dence on the coupling constants. By comparing these 
results with the dispersion curves for SrCuC>2 reported 
in Ref. 3 we can estimate the Hubbard effective coupling 
constants appropriate of this material: t ~ 0.53 eV and 
U/t ~ 7. 



III. THE QUASI-PARTICLE WEIGHT FROM 
LANCZOS DIAGONALIZATION 



Unfortunately, the formal Bethe Ansatz solution does 
not lead to a practical way for the evaluation of the quasi- 
particle weight (5) at arbitrary couplings and therefore 
we resort to Lanczos diagonalizations in lattices up to 
L = 14 sites. As previously noticed, we first have to de- 
vise a method to select the correct single spinon states at 
finite coupling U/t, for these states are not identified by a 
good quantum number at finite U. Our method is based 
on an adiabatic procedure starting from the strong cou- 
pling limit. We first perform a Lanczos diagonalization 
on the (L — l)-sitc Hciscnberg model in the symmetry 
subspace of total momentum Q, leading to the numerical 
determination of <pH and of the exact eigenstates of the 
one-hole Hubbard model for U — > oo via Eq. (3). In this 
limit, the single spinon states are indeed the lowest en- 
ergy eigenstates at fixed spinon momentum Q and can be 
easily obtained by Lanczos (or DMRG) technique, whilst 
at finite U the relevant intermediate states are not nec- 
essarily in the low excitation energy portion of the Hub- 
bard spectrum. Then we take advantage of the continu- 
ity of the one spinon states between the weak and strong 
coupling limit by adiabatically lowering the interaction 
strength U and performing successive Lanczos diagonal- 
izations for smaller and smaller couplings U n . At the n th 
step we keep the exact eigenstate having the largest over- 
lap with the eigenstate at the (n— l) th level. In this way 
we are able to identify the single spinon states down to 
small values of U ~ t, each state being uniquely identified 
by Q, i.e. by the momentum of the "parent" Hciscnberg 
eigenstate. 

A check on the validity of the singlc-spinon approx- 
imation comes from the completeness condition on the 
intermediate states: 

ni(k) = (*ol4 iiCM |* ) = E K*iKil*o)| 2 

{|*i>} 

> £|(-fc,Q| CM |vI/ )| 2 



(7) 



where n^{k) is the momentum distribution of the down 
spins at half filling and the equality holds if and only if 
the single spinon states included in the sum via the defi- 
nition of the quasi-particle weight Zk(Q) (5) exhaust the 
spectral weight at each k. The amount of violation of 
this sum rule quantifies the weight of all the neglected 
states in the Hilbert space due to the single-spinon ap- 
proximation. In Fig. 3 we plot nj.(fc) and TIq 2>k(Q) 
restricted to the one spinon states: the violation of the 
completeness condition is smaller than 0.01 at all fc's 24 . 
Note how, even at fairly large values of U/t, the momen- 
tum distribution is considerably depressed for k larger 
than the free electron Fermi momentum kp = ^, strongly 
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FIG. 3: Momentum distribution rij(fc) (black open circles) 
and -j^-j- YIq Zk(Q) (red triangles) versus the hole momentum 
k from Lanczos diagonalization, for U = It in a L — 14 ring. 




FIG. 4: Panel (a): Zk(Q) versus spinon momentum Q for 
(7 = 1004 and different lattice sizes (x: L=6, A: L=8, ■: 
L=10, *: L=12, •: L=14). Open (green) symbols for total 
momentum k — n/4 and full (black) symbols for k = 0. Lines 
are polynomial fit to Lanczos data. Skewed boundary con- 
ditions are used in order to fix the same total momentum of 
the state — k (relative to half filling) for all L's. Panel (b): 
same as (a) for U = It. Panel (c): binding energy at fixed 
k referred to half filling versus spinon momentum Q in the 
thermodynamic limit E = Eo — E\ for U = 1004. Panel (d): 
same as (c) for U = 7t. 



reducing the spectral weight in the second half of the Bril- 
louin zone. This feature is consistent with the photoemis- 
sion experiments performed with high energy photons 3,6 . 
Conversely, in the strong coupling limit U — * oo, the mo- 
mentum distribution becomes flat, n±(k) = 1/2, washing 
out this effect. 

The dependence of the quasi-particle weight on the 
strength of the Coulomb repulsion has been investigated 
and is summarized in Fig. 4 for strong and intermedi- 
ate U/t and for different lattice sizes. The quasi-particle 
weight has been evaluated by Lanczos diagonalization on 
lattices ranging from L = 6 to L = 14 sites. By using 
standard periodic boundary conditions, the total momen- 
tum of the state would be quantized in units of 2tt/L, 
making size scaling impractical. In order to avoid this 
problem we have adopted skewed boundary conditions: 



given an arbitrary hole momentum k we choose the flux 
at the boundary in such a way to match k with the quan- 
tization rule. Figure 4 reveals an astoundingly negligible 
size dependence and the expected vanishing of the quasi- 
particle spectral weight outside the spinon Fermi surface, 
with singularities at the Fermi momenta. While Zk(Q) 
is almost independent on k at large U, as expected 8 , it 
shows more structure for realistic values of U ft. The fur- 
ther peak (or shoulder) present for k < \ is indeed remi- 
niscent of the free Fermi nature of the electrons at U = 0. 
In the free particle limit, only one state provides a finite 
contribution to the spectral function: the holon sits at 
the bottom of the band (fc/, = 7r) and the quasi-particle 
weight Zk(Q) reduces to a delta function at Q = n + k. 
When such a form of Zk(Q) is substituted in Eq. (4), 
the known free particle result is recovered. Remarkably, 
a remnant of the free particle peak in Zk (Q) is still visible 
at U = 7t, as shown in Fig. 4b. 

IV. WEAK COUPLING LIMIT 

The Green's function of one dimensional models has 
been thoroughly investigated by bosonization methods: 
while in the Luttinger liquid regime its asymptotic form 
is characterized by power law tails 1 , precisely at half fill- 
ing the Green's function is known to display a more com- 
plex behavior due to the presence of a gap in the holon 
spectrum. At weak coupling, the holon dispersion near 
the bottom of the band shows a "relativistic " structure: 

e h (k h ) = ^v 2 h 6k 2 h +m 2 (8) 

where m is the charge gap and 8kh = kh ± tt is the holon 
momentum measured from the bottom of the band. Note 
that the holon spectrum (8) is shifted by fjb = U/2 with 
respect to our previous definition. In Fig. 5a we plot 
the exact Bethe Ansatz spectrum at U = 3t and the 
form (8) predicted by bosonization with suitably chosen 
parameters m,Vh- 

In order to compare the results of our single spinon ap- 
proximation with the bosonization form, it is convenient 
to introduce the single hole Green's function in imaginary 
time: 

G;(fc, r) =< tt |4 4 Cktl \* > (9) 

According to bosonization 9 , the Green's function 
Gi(k,T) of a hole of momentum close to kp = ir/2 ac- 
quires a factorized form in real space: 

Gf{x,r) ^ J ^Gf(k,r)e^ 

= e^ x G h (x,T)G s {x,T) (10) 

where the superscript R identifies the contribution to the 
Green's function due to right moving holes. Here, Gh and 
G s just depend on holon and spinon degrees of freedom 
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respectively. The spinon term is simply given by 
G s (x,t) 



y/v s T + IX 



(11) 



while the holon contribution is predicted, by the form 
factor approach, to behave as 

/ \ Fni f \q/2 — mr cosh 6 — im — sinh d\ 

Gh(x,T)=T\ — / d8ei J 

V Vh J-oo 

(12) 

with T - 0.0585. .. 9 . The question now arises whether 
our single spinon approximation is consistent with such 
a factorized form. By inserting a complete set of inter- 
mediate states into the definition (9) and adopting the 
single spinon approximation, the full Green's function in 
imaginary time can be written as 



Q 



,-[«h(feh)+e,(Q)]i 



(13) 



where the momentum conservation relation fc = kh + 
Q is understood and the holon spectrum e^(fc/ l ) is now 
referred to the chemical potential fi. Notice that, due to 
momentum conservation, the combined requirements of 
having k ~ kp = ^ and kh ~ — 7r (i.e. the hole sits near 
the bottom of the band) force Q ~ By substituting 



the asymptotic forms (8) and e s 
we get: 



'37T 

• 2 



1 



v s q for q > 



G?{x,t) = e^ x — [* **e-to**"ri 



dO 

—Z k (Q) cosh e 

Z7T 



(14) 

— im — sinh 6 — rar cosh 6 



where we set Skh = ——sinh 9. This form does indeed 

Vh 

factorize in a holon and spinon part, as predicted by 
bosonization, provided the quasi-particle weight does: 



Z k (Q) ~ Zh{k h )Z s {Q) 



(15) 



Notice that our approach, being based on a numerical 
evaluation of the quasi-particle weight, docs not allow for 
an independent demonstration of such a factorized form. 
We just observe that the bosonization approach and the 
single-spinon approximation lead to the same result if we 
assume that Eq. (15) holds. Following Ref. 14 we argue 
that the factorization of the quasi-particle weight at low 
energies (15) reflects the trivial structure of the holon- 
spinon scattering matrix in this limit. 

As previously noticed, the spinon contribution to the 
quasi-particle weight gives rise to the square root diver- 
gence at the spinon Fermi surface, with leading behavior 
Z s (4^ — q) ~ q^ 1 / 2 for q > which correctly reproduces 
Eq. (11) when the ultraviolet cut-off in Eq. (14) is dis- 
regarded. Matching Eqs. (12) and (14) then selects a 
unique form of the holon quasi-particle weight: 



Zh{k h ) 



\/th{kh) - VhSk h 
th(kh) 



(16) 
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FIG. 5: Panel (a): Holon spectrum eu{kh) of the Hubbard 
model at U = 3t from Bethe Ansatz (shifted by /i = U/2) 
compared to the Lorentz form (8). The fitting parameters are 
m = 0.316t and Vh = 2.43i. Panel (b): Dimensionless holon 
quasi-particle weight Zh = . —Z^ in the single spinon ap- 



proximation for U — 3t and lattice sizes L — 10 (full squares) , 
L — 12 (empty circles), L — 14 (full circles), compared to the 
bosonization result (16) (line). 



with eh(kh) given by Eq. (8). The scale factor in Eq. 
(16) has been fixed by evaluating the Green's function in 
the r — > limit, where it coincides with the momentum 
distribution. In Fig. 5b we compare Eq. (16) with the 
numerical results for Zhikh) obtained by Lanczos diag- 
onalization at U = 3t. No fitting parameters have been 
used: In order to obtain Zh(kh) we first evaluated Z^{Q), 
as discussed in Section III, then we divided the result by 
Z S {Q) ~ q~ x l 2 evaluated at the spinon momentum Q 
closest to the spinon Fermi point Qp = The two 
parameters Vh and m arc independently obtained from 
the holon spectrum (also shown in Fig. 5a). As usual 
the Lanczos data display a very small size dependence 
and allow for a precise identification of the holon quasi- 
particle weight Zhikh)- The agreement between the two 
expressions is remarkable for Skh > while some discrep- 
ancy is found for negative Skh- Note however that the 
asymptotic form of the holon quasi-particle weight (16) 
holds only at low energies and weak coupling, while the 
comparison shown in Fig. 5 is performed for U = 3t. 
The results at lower values of U/t are plagued by severe 
finite size effects: in the U — ► limit, the holon mass 
m vanishes exponentially and the dimensionless momen- 
tum scale m/vh vanishes as well. Therefore, at very weak 
coupling, the relevant holon momenta are constrained in 
an extremely small interval around kh = ±tt, a range not 
easily accessible due to the momentum quantization rule 
in finite Hubbard rings. 



V. RESULTS FOR SrCu0 2 

We are now ready to compare our results for the spec- 
tral function of the ID Hubbard model with precise pho- 
tocmission data recently obtained for SrCu02 3 . A pre- 
liminary study, based on the strong coupling limit of 
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FIG. 6: Locus of singularities of the spectral function for 
the Hubbard model at U — 3.7 eV and t — 0.53 eV (open 
circles) compared to the experimental results by Kim et al. s 
(full circles) for k± = 0.1. 



the Hubbard model, pointed out some discrepancies, re- 
lated to the peak heights and widths 3 . Fig. 6 shows 
the singularity loci of the ID Hubbard model with the 
parameters t = 0.53 eV and U = 3.7 cV, together 
with the ARPES results from Kim et al 3 . The nice 
agreement suggests that this material indeed represents 
a good experimental realization of the simple one dimen- 
sional Hubbard model. The effects due to inter-chain 
coupling, phonons, finite temperature and other pertur- 
bations appears rather small and mostly limited to the 
spinon branch. We remark that the same material has 
been already theoretically investigated on the basis of the 
Hubbard and t-J model by several groups 3 ' 25,26 leading 
to different sets of parameters both for the hopping in- 
tegral 0.3eV < t < 0.7eV and for the Coulomb repulsion 
2eV < U < 6.5eV. Our analysis shows that both spin 
and charge fluctuations play a key role in determining 
the line shape of the spectral function of the Hubbard 
model, even at moderately high values of the coupling 
U/t. 

In Fig. 7 the spectral function has been plotted versus 
the binding energy E = —lu for three representative val- 
ues of the total hole momentum k. The experimental line 
broadening reported in Ref. 3 has been also included in 
the Hubbard model results, leading to a merging of close 
peaks. The density plot clearly reproduces the overall 
shape defined by the singularities of the spectral function 
shown in Fig. 6. As expected, most of the spectral weight 
is indeed concentrated in the first half of the Brillouin 
zone between the holon and the spinon band. Although 
the relative intensity of the ARPES signal at the two sin- 
gularities depends on the details of the band structure, 
the power-law nature of the divergences implies that the 
intrinsic width of each peak is always comparable with 
the separation between the holon and the spinon branch 
Auj ~ eh(k + 77) — [e/j(7r) + 6 s (tt + k)]. The average in- 
tensity can be estimated on the basis of the sum rule 
(7) and scales as (Aw)" 1 ' 2 , getting smaller when the two 
branches separate, as shown both in experiments 3 and in 
numerical calculations 7 . 
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FIG. 7: Upper panels: A(k,u) calculated via Eq. (4) for three 
representative values of the momentum k and U = It. The 
binding energy is E — —u). The lower curves are the convo- 
lution of the spectral function with a Gaussian with FWHM 
equal to 0.12i corresponding to an experimental resolution of 
60 meV 3 . Lower panel: density plot of the spectral function 
on the (k, E) plane with the experimental broadening. 



VI. CONCLUSIONS 

The single spinon approximation, combined to Bcthe 
Ansatz results and Lanczos diagonalizations allows to 
obtain very accurate results for the dynamic properties 
of a single hole in the one dimensional Hubbard model. 
The Lehmann representation of the spectral function (2) 
shows that two separate ingredients combine to define the 
overall shape of A(k,ui): the excitation spectrum and 
the quasi-particle weight. The idea at the basis of our 
method is to limit the size effects that plague numerical 
results by dealing with these two quantities separately: in 
the ID Hubbard model the excitation spectrum is given 
exactly by the Bethe Ansatz equations in the thermody- 
namic limit 11 , while the quasi-particle weight is obtained, 
in the single spinon approximation, by Lanczos diagonal- 
ization. Size effects are shown to be negligible and the 
accuracy of the approximation can be checked a posteri- 
ori by a frequency sum rule (7). Our expression for the 
spectral function of the ID Hubbard model (4) is con- 
sistent with the structure predicted by bosonization at 
weak coupling, provided the quasi-particle weight Zk(Q) 
factorizes as shown in Eq. (15). A numerical test carried 
out at U = M docs not show a convincing quantitative 
agreement with the result obtained by the form factor 
approach 9 , possibly due to the difficulty to achieve the 
U -> limit. 

The extension to finite doping is straightforward but 
in principle this approach can be also generalized to 
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other fermionic lattice models, in one or more dimen- 
sions, provided the relevant states entering the quasi- 
particlc weight in the strong coupling limit can be easily 
classified. This would be the case of the extended Hub- 
bard model (i.e. a Hubbard model with nearest neighbor 
Coulomb repulsion) or in the presence of lattice dimer- 
ization. Clearly, for non integrable models, no analytical 
information on the excitation spectrum is available and 
a size scaling on the energy spectrum is also required. A 
study of such generalizations may be useful to understand 
the role of some perturbation on the spectral function of 
correlated electron models. 

The specific example of SrCU02 shows that our 
method allows for a direct comparison between theory 



and ARPES experiments and for an accurate determi- 
nation of the Hubbard parameters which best describe 
the hole dynamics in the material. The spectral func- 
tion derived here provides a natural explanation of the 
observed reduction of the spectral weight in a half of the 
Brillouin zone and of the broad line-shape detected in ex- 
periments. Future applications of this method to the case 
of cold atoms in optical traps may help in pointing out 
the peculiar features of one dimensional physics in other 
experimental realizations of correlated one dimensional 
Fermi gases. 

We thank C. Kim and F.H.L Essler for stimulating 
correspondence. 
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